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ABSTRACT 


We present a precise measurement of downward-going albedo proton fluxes 
for kinetic energy above ~ 70 MeV performed by the PAMELA experiment at 
an altitude between 350 and 610 km. On the basis of a trajectory tracing simula¬ 
tion, the analyzed protons were classihed into quasi-trapped, concentrating in the 
magnetic equatorial region, and un-trapped spreading over all latitudes, including 
both short-lived (precipitating) and long-lived (pseudo-trapped) components. In 
addition, features of the penumbra region around the geomagnetic cutoff were in¬ 
vestigated in detail. PAMELA results signihcantly improve the characterization 
of the high energy albedo proton populations at low Earth orbits. 


1. Introduction 


The Cosmic Ray (CR ) albedo radiation has been detected since late 1940 s by rockets 


and balloon experiments fjVan Allen &: Ganges 


1956 : 


Ormes fc Webber 


1964 ). These 


particles, generated with upward-going directions by the interaction of CRs from the 
interplanetary space with the Earth’s atmosphere, are typically classihed into “splash” and 
“re-entrant” alb edo componen ts, depending on the type of trajectories followed from the 


production site flTreiman 


19531) : while the former refers to particles which are able to escape 


from the magnetosphere (i.e. allowed Stormer trajectories), the latter includes particles 
whose trajectories are bent by the geomagnetic held back to the Earth (i.e. forbidden 
trajectories). 

The re-entrant albedo (hereafter albedo) population comprises quasi-trapped and 
un-trapped components. The fo rmer concent r ates in the near equatorial region, inside ari d 


below the inner Van Allen belt flMoritz 


1972 


Hovestadt et al. 


19721: 


Fiandrini et al. 


29941) ■ 
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and presents featnres similar to those of stably-trapped particles in the radiation belts 


althongh it is characterize d by limited li 


spreads ove r all latitndes 


penumbra fjStormer 


1950; 


Alcaraz et ah 


Shea et al 


etimes and much 


2000 


Bidoli et ah 


ess intense fluxes. The latter 


20031) including the so-called 


19651) ■ a region around the local geomagnetic cutoff 


rigidity where particle trajectories be come chaotic and both allowed and forbidden bands 


of CR particle access are present (see 


Cooke et ah 


1991 


for dehnitions of CR cutoffs). 


orbits have been reported by the PAMELA mission flAdriani et al 


> 70 MeV) Cl 

h radiation at low Earth 

Adriani et al. 

2014 

). Thanks to the 


orbit and the high identihcation capabilities, the instrument is able to provide detailed 
information about particle fluxes in different regions of the terrestrial magnetosphere, 
including energy spectra, spatial and angular distributions. PAMELA measurement of 
geomagnetically trapped protons, and qua si-trapped electrons and positrons can be found 


in publications flAdriani et al. 


2015 


2 OO 9 I) . In addition, for the hrst time P AMELA has 


revea. 


ed the existence of a trapped antiproton component in the inner belt flAdriani et al 


2nilal ). In this article the measurement of the albedo proton fluxes is presented. 


2. PAMELA Data Analysis 


PAMELA is a space-based experiment designed for a precise measurement of the 


charged cosmic ra diation in the 


hundreds of GeV fjPicozza et al. 


dneti c energy range from some tens of MeV up to several 

brmance, particle 


20071) . Details about apparatus per: 


selec 


ion and selection efficiencies evaluation can be found elsewhere flAdriani et ah 


2011b 


20131 ). The Resurs-DKl satellite, which hosts the apparatus, has a semi-polar (70 deg 


inclination) and elliptical (350 610 km altitude) orbit, and it is three-axis stabilized. The 

spacecraft orientation is calculated by an on board processor with an accuracy better than 
1 deg which, together with the good angular resolution (< 2 deg) of the PAMELA tracking 
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system, allows particle direction to be measured with high precision. The PAMELA major 
axis is mostly oriented toward the zenith direction. 


2.1. Data Set 


The analyzed data set includes protons acquired by PAMELA between 2006 July and 
2009 September, corresponding to a total live time of about 800 days. Data registered 
between 2006 December 6 and 2007 January 10, characterized by high geomagnetic 
isturbance levels as sociated with the December 6, 13 and 14 solar particle events 


flAdriani et ah 


2011a), were excluded from the sample. The time variation of the PAMELA 


detector performance, mostly due to the sudden failure of some front-end chips in the 
tracking system, were carefully studied and accounted for in the efficiencies estimate. 


2.2. Geomagnetic Models 


The IGRFIO ( Mans fc Macmillanll2010 1 and the TS05 f Tsvganenko fc SitnovI 12005 1 


models were used for the description of internal and external geomagnetic held sources, 
respectively: the former is a spherical harmonic expansion of the Earth’s main magnetic 
held; the latter is able to provide a dynamical description (with a hve-minute resolution) of 
the geomagnetic held in the inner magnetosphere, based on recent satellite measurements. 


Solar wind and inter planetary magnetic 


Omniweb database flKing fc Papitashvili 


leld v alues were derived from the high resolution 


20041 ). All parameters were evaluated on an 


event-by-event basis through linear interpolation. 
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2.3. Coordinate Systems 


Data were analyzed in terms of “Altitude Adjusted Corrected GeoMagnetic” (AACGM) 
coordinates, developed to provide a more realistic description of high latitude regions by 
accounting for the multipolar geomagnetic held. They are dehned such that all points 
along a magnetic held line have the same geomagneti c latitude and longitude, so that they 


are c l osely related to invari ant magnetic coordinates flBaker et ah 


1989 


Gustafsson et ah 


1992 


Heres fc Bonito 


2007). The AAGGM reference 


“Gorrected GeoMagnetic” (GGM) coordinate system 
surface. 


rame coincides wi 


Brekke et al. 


h the standard 


1997 ) at the Earth’s 


3. Trajectory Reconstruction 


Using spacecraft ephemeris data (position, orientation, time), and the particle rigidity 
[R = momentum/charge) and direction provided by the magnetic spectrometer, trajectories 


of all selected protons were reconstructed in the Earth’s magnetosphere 


tracing program based on numerical integration methods (jSmart fc Shea 


3Y means o 


2000 


20051) . 


implementing the aforementioned geomagnetic held models, and carefully optimized for the 
calculation of very long trajectories in the magnetosphere. 


Measured rigidities were corrected for the mean energy loss in the apparatus, evaluated 
with realistic simulations of the instrument. This solution was preferred to standard 
statistical methods (e.g. huxes unfolding) in order to provide a more reliable input to the 
trajectory simulation, and to preserve spectral structures especially near the penumbra 
region. The ehects of the hnite rigidity and angular resolution of PAMELA on the trajectory 
reconstruction were accounted for in the estimate of systematic uncertainties. 




























7 


3.1. Proton Classification 


Trajectories were propagated back and forth from the measurement location, and 
traced until: 


1 . they reached the model magnetosphere boundaries; 

2 . or they intersected the absorbing atmosphere limit, which was assumed at an altitude 
of 40 km (mean albedo production altitude); 

3. or they performed more than ~ steps for both propagation directions, where 

R is the particle rigidity in GV. 


Particles satisfying the hrst condition were classihed as protons from interplanetary space. 
The second case corresponds to albedo protons. Finally, the third category includes 
stably-trapped protons from the inner Van Allen belt, detected in the so-called “South 
Atlantic Anomaly” (SAA) region at PAMELA altitudes. Since the program uses a dynamic 
variable step length, which is of the order of 1% of a particle gyro-distance in the magnetic 
field, the applied rigidity-dependent criterion (3) ensures that at least four drift cycles 


around the Earth were 

)erformed: the measurement of such a component is described 

elsewhere ( 

Adrian! et al. 

2015 

). The initial sample also included a negligible fraction of 


escaping albedo protons with rigidity above the geomagnetic cutoff. 


Each albedo trajectory was checked by considering, as a first approximation, the 
particle motion in the magnetosphere as a superposition of three quasi-periodic motions: a 
gyration around field lines, a bounce between magnetic mirrors in opposite hemispheres, 
and an azimuthal drift around the Earth. By using combined selections on corresponding 
mean frequencies {ugyro, ojhounce and Udrift) as a function of several variables of interest 
(energy, pitch angle, latitude, etc.), the selected sample was sub-divided into two main 
categories: 





• events with trajectories similar to those of stably-trapped protons, bnt originated 
and re-absorbed by the atmosphere dnring a time shorter than a few drift periods, 
were identihed as qnasi-trapped (QT). Their trajectories were verihed to satisfy the 
adiabatic conditions, in particular the hierarchy of temporal scales: ojhounce/^gyro 

< ei and ujdrift/^^bounce ^ ^ 2 , where ei and 62 were estimated to be of the order of 

10-2 ^ 10-h 

• The rest of the sample was classified as un-trapped (UT). Qualitatively, two 
sub-components can be identified: 

— precipitating protons (UTs), with lifetimes shorter than a bounce period. 
Corresponding Ubounce values are similar to those of quasi-trapped protons, while 
ojgyro distribution is much broader outside the SAA, extending to much lower 
values. 

— Pseudo-trapped protons (UT^,), with relatively long lifetimes. They are 
characterized by large gyro-radii and udrift, and by small Ugyro and Ubounce 
values, resulting in unstable trajectories due to resonances occurring between 
component frequencies. They can perform several drift cycles (up to a few 
hundreds) reaching large distances from the Earth’s surface, sometimes forming 
intermediate loops, before they are re-absorbed by the atmosphere. 


3.2. Lifetime Distributions 

The lifetime of the different populations was estimated from the tracing simulation 
as the time between the particle origin (traced backward) and its subsequent absorption 
(traced forward) in the atmosphere (i.e. the tracing time r). Results are displayed in 
Figure dl where the lifetime for the measured sample is shown as a function of the AACGM 


9 


latitude A (top-left panel), and of the product R ■ L 


(top-right pan el), where R is the 


19661) in units of Earth’s 


particle rigidity and L is the Mcllwain’s parameter flMcIlwainl 
radii (Rk); note that R ■ 14.3 GV Rg corresponds to the Stormer vertical cutoff 


(IStormer 


1950 


Shea et ah 


19651 ) for the PAMELA epoch. Distributions as a function of 


kinetic energy are also reported in bottom panels, for magnetic latitudes |A| <35 deg (left) 
and |A| >35 deg (right) respectively. 

The several albedo components can be easily distinguished. QT proton lifetimes 
typically amount to r ~ 0.3 30 s, while the upper bands (r up to ~ 2 min) correspond to 

particles performing more than a revolution around the Earth. Si nce the QT proton lifetime 


is of the order of a half drift period, which scales with 1 / 7 / 5 ^ flWaltlll994l ). lifetime and 
energy have an approximately inverse proportional relation. 


Conversely, the UT 5 proton lifetime ranges from some fraction of s up to a few s 
(depending on latitude) and it is shorter than the typical bounce period, which scales with 
1//3, resulting in a weaker dependency on energy. In particular, their distribution is given 
by two superimposed bands, corresponding to particles crossing the magnetic equator once 
and twice, respectively. In addition, a population of very-short-lived (< 30 ms) protons 
emerges at high energies (0.8 8 GeV) and low latitudes (|A| < 10 deg), with a tracing time 
increasing with energy: because of the large gyro-radius with respect to the magnetic held 
curvature, such particles have no stable bounce motion {ubounce ~ ^gyro, large udrift values). 

Non-adiabatic or large gyro-radius effects cause the breakdown of (quasi-) trapping 
conditions for i? ■ >4 GV R|.. Such a geomagnetic region is populated by UT l 

protons (see top-right panel), characterized by irregular trajectories with no periodicity; 
consequently, their lifetime distributions are much broader, extending up to a few tens of 
min (pseudo-trapped motion). Such a component can be energetic enough to populate the 
penumbra region (i? ■ > 10 GV R|;), where particles of both atmospheric and galactic 
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origin are present. 


3.3. Production and Absorption Points on the Atmosphere 


The production and absorption points on the atmosphere are shown in Figure O 
respectively in left and right panels. In the case of QT protons (upper panels), the former 
are located in a region extending westward from the SAA, where the Earth’s magnetic held 
is at its minimum due to asymmetries between terrestrial magnetic dipole and rotational 
axes (tilt ~10 deg, offset ~500 km). While drifting from the SAA, protons encounter 
stronger magnetic helds and the altitude of their mirror points increases, until they reach 
again weaker magnetic held regions; then their mirroring altitude decreases and hnally 
they are absorbed by the atmosphere, mainly on the region on the East side of the SAA. 
Both source and sink points are located in two regions, in the southern and in the northern 


magnetic hemi sphere respec 


magnetic held flHuang et ah 


dvely, as a consequence of the multipole moment of the Earth’s 


200l|). 


Conversely, since UT 5 protons (middle panels) are created and absorbed by the 
atmosphere in a very short time, their production and absorption points are located 
near the detection point, uniformly populating the whole geomagnetic region explored by 
PAMELA. In particular, UT^ sink points have a peak in the SAA, while origin points have 
two peaks: the former in the SAA, the latter in the northern magnetic region corresponding 
to southern mirror points in the SAA. Particles produced in the SAA can undergo a 
rehection in the upper hemisphere (where mirror point altitude is higher) before being 
re-absorbed, and correspond to the upper band (r larger than a half bounce period) in the 
UT 5 lifetime distributions (see Figured]). 


Similar to UT 5 protons, distributions (bottom panels) spread over all longitudes. 
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In addition, their absorption points have a minimum in the SAA, while source distribution 
has a minimum in the geomagnetic region opposite to the SAA (around ~20° N, ~90° E). 
It should be noted that reported distributions were evaluated at a hxed altitude of 40 km, 
used as a limit for the absorbing atmosphere, while the actual source and sink positions are 
expected to be broader. 


4. Albedo Proton Fluxes 

Fluxes were derived for bins of AACGM longitude, latitude and proton energy, by 
assuming an approximately isotropic particle distribution. Integral flux maps (m“^ s“^ sr“^) 
in geomagnetic coordinates are displayed in Figure |3] for different energy bins. Results for 
the different populations are reported. The flux dependency on longitude is quite complex, 
reflecting the asymmetries between the magnetic and geographic axes, and it is signihcantly 
influenced by the presence of non-dipolar components of the geomagnetic held. In addition, 
at PAMELA energies the particle gyro-radius is quite large, affecting the portion of huxes 
which can be investigated by PAMELA, so that measured distributions have different 
shapes for different energies. Low energy QT and UTg distributions display a peak near 
the SAA (longitude ~ 0 50 deg), where the Earth’s magnetic held is at its minimum. In 

addition, we found that measured energy spectra become harder while approaching to the 
geomagnetic region opposite to the SAA, where the geomagnetic held has a local maximum. 
Fluxes were averaged over the pitch angle and the altitude ranges covered by PAMELA. An 
average hux increase of a factor ~ 1.5 was observed between lowest and highest explored 
altitudes. 
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4.1. Energy Spectra 


Energy spectra were evaluated inside and outside the SAA, by requiring local magnetic 
field intensities lower and higher than 0.23 G, respectively. Figure 0] shows the spectra 
outside the SAA measured at different AACGM latitudes |A|. Fluxes were averaged over 
longitudes. Reported error bars account for statistical and systematic uncertainties, and 
are smaller than the plotted points in most cases. The several proton components (QT, 
UT 5 and UT^,) are denoted with different colors (blue, green and red, respectively); galactic 
proton spectra (black points) are also reported for a comparison; finally, the underlying 
orange points denote the total (albedo+galactic) proton spectra. 

Fluxes of QT protons are limited to low latitudes and to energies below ~ 8 GeV; 
they smoothly decrease with increasing latitude and energy. Gonversely, UT 5 distribution 
spreads to higher latitudes, with energetic spectra extending up to ~10 GeV. Their flux is 
suppressed at the equator, since corresponding magnetic field lines typically do not reach 
PAMELA altitudes. Both QT and UT 5 components have a low energy peak in the SAA, 
while spectra measured outside the SAA are harder. 


Finally, the UT^ protons concentrate at the highest latitudes and energies (up to ~ 20 
GeV), with a peak in the penumbra. Such a feature is due to large gyro-radius (10^ 10^ 

km) effects: particles with increasing rigidity come from a more and more extended region 
causing fluxes to increase, until trajectory curvature becomes tog larg e and they escape 


from the magnetosphere, and fluxes decrease again flBobik et al 


20061) . Gorresponding 


pitch angle distributions are much broader, extending to highest values. Data supporting 
Figure m are available in the Supporting Information Tables S1-S4. 


Figure [5] reports the comparison between total proton spectra outside the SAA {B > 
0.26 G), including both galactic and albedo components, measured by PAMELA (red) and 
AMS-01 (black), for different bins of geomagnetic latitude |A|. Note that AMS fluxes were 
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evaluated (using CGM coordinates flBrekke et al 


1997)) at altitudes of 350 390 km and 


are limited to geographic latitudes lower than ~ 52 deg. Latitudes are reported in radians, 
for consistency with AMS published results. Despite good qualitative agreement, signihcant 
differences can be observed. AMS spectra results up to 50% lower than PAMELA ones, 
especially at low energies and latitudes. On the other hand, AMS fluxes are a factor 2 
3 larger than PAMELA ones near the penumbra. In addition, PAMELA spectra show 
some structures near the dip region before the penumbral peak, especially at middle-low 
latitudes, which are absent in AMS data. 


4.2. Penumbra Protons 


Features of the penumbra region can be more deeply investigated in Figure [6l where 
the fraction of albedo protons is displayed as a function of particle rigidity and AACGM 
latitude (left panels); for a comparison, distributions as a function of Mcllwain’s L-shell 
are also shown (right panels). Bottom panels report corresponding rigidity projections, 
including statistical uncertainties on measured distributions. The penumbra was identified 
as the region where PAMELA detected both albedo and galactic proton trajectories. 

The black curves denote a ht of points with an equal percentage of the two components: 


77(*%(A) ~ 11.62 • cos^A - 0.215 GV 

and 

~ 14.14/L2 - 0.187 GV 


for magnetic latitude A and L-shell distributions, respectively. For a comparison, 
the red curve indicates the Stormer vertical cutoff evaluated for the PAMELA epoch: 
R<^v(L) ~ 14 . 3 GV. The htting functions 7?5o% are similar to the one used by 


flOdiore et al. 


200ll) to reproduce cutoff observations made by the MAST instrument on 
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the SAMPEX satellite, and they account for magnetospheric effect s at high latitudes. 


However, because of the non-dipole terms in the geomagnetic held fjSmart &: Shea 


200511 ■ 


the ht results can vary up to a factor 10 % at equatorial and mid-latitude locations. 

The usage of a realistic external geomagnetic held model has a signihcant impact on 
this study inhuencing the reconstructed fraction of albedo protons, as highlighted in Figure 
[3 The top panel compares the estimated values derived from PAMELA data 

with and without the inclusion of external held sources, corresponding to the red (IGRFIO 
-|- TS05) and the blue curve (IGRFIO) respectively. In the latter case, we found that ht 


results {RlliL) - 14.22/L^ GV 


are in a very good agreement with the standard cutoh 


relation Rsv{L) fjShea et al.lll985l ) at all L-shells. On the other hand, signihcant diherences 


(up to a factor ~ 2) can be noted between the IGRFIO model and the combined model 
conhguration at high L, as shown in the bottom panel. Similar results were obtained for 
the A distribution. 


4.3. Albedo Fluxes in the SAA 


Figure [H] shows the average energy spectra of QT, UT 5 and UT^, protons in the SAA. 
Measured spectra are softer than those in other magnetic regions, with a signihcant hux 
increase at low energies. In case of QT and UT 5 protons, an additional contribution is 


possibly expected from the decay of free albedo neutro ns, whic 


1 constitutes the main source 


of geomagnetically trapped protons in radiation belts (jSinger 


1958 


Farley fc Walt 


197l|). 


The UTi, component emerges at high energies (> 1.5 GeV), concentrating near the local 
geomagnetic cutoh. 
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5. Conclusions 

PAMELA measurements of energetic (> 70 MeV) albedo proton fluxes at low Earth 
orbits (350 -i- 610 km) have been presented. The detected sample, corresponding to 
data acquired by PAMELA between 2006 July and 2009 September, was analyzed and 
classihed on the basis of a realistic simulation of particle trajectories in the Earth’s 
magnetosphere into quasi-trapped and un-trapped protons: the former consists of relatively 
long-lived protons, detected in the near equatorial region, with trajectories similar to 
those of stably-trapped protons from the inner radiation belt; the latter was found to 
spread over all explored latitudes, including short-lived (precipitating) protons together 
with a long-lived (pseudo-trapped) component constituted by particles with rigidities 
near the local geomagnetic cutoff and characterized by a chaotic motion (non-adiabatic 
trajectories). Eluxes were mapped by using the AACGM coordinates to provide a more 
realistic description at higher latitudes. 

PAMELA results signihcantly enhance the characterization of high energy albedo 
proton populations in a wide geomagnetic region, including the SAA and the penumbra, 
enabling a more accurate and complete view of atmospheric and geomagnetic effects on the 
CR transport near the Earth. 
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Fig. 1.— Top panels: lifetime distributions for the measured sample (counts) as a function 
of the AACGM latitude A (left) and of the product of particle rigidity R and L-shell squared 
(right). Bottom panels: distributions as a function of kinetic energy, for AACGM latitudes 
|A| lower and greater than 35 deg (left and right, respectively). The several albedo com¬ 
ponents can be easily distinguished: quasi-trapped (QT), precipitating (UT 5 ) and pseudo- 
trapped (UT^,). 
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Fig. 2.— Count distributions of the production (left panels) and absorption (right panels) 
points on the atmosphere (40 km) for quasi-trapped (top), precipitating (middle) and pseudo- 
trapped (bottom) protons as a function of the geographic coordinates. 
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Fig. 3.— Proton integral fluxes (ni“^ sr“^) as a function of AACGM longitude and 

latitude, for different energy bins. Results for the several albedo populations are reported 
(from left to right): quasi-trapped (QT), precipitating (UT 5 ) and pseudo-trapped (UT^) 
protons, along with the total sample. Note that a different energy binning is used for QT 
protons (left column). 
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Fig. 4.— Differential energy spectra (GeV“^ m~^ sr“^) outside the SAA region mea¬ 

sured for different bins of AACGM latitude |A| (see the labels). Results for the several 
proton populations are shown: quasi-trapped (QT, blue circle), precipitating (UT^, green 
squares), pseudo-trapped (UT^, red triangles) and galactic (black empty circles); the under¬ 
lying orange points denote the total proton flux. Lines are to guide the eye. 
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Fig. 5.— Total proton spectra, including both galactic and albedo com ponents, for different 
bins of geomagnetic latitude |A| (radians). For comparison, AMS-01 flAlcaraz et al.ll2000[) 
results (black points, |A| < 1.1 rad) are also reported. Lines are to guide the eye. 
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Fig. 6.— Top panels: fraction of albedo protons in the pennmbra region, as a function of 
particle rigidity and AACGM latitnde |A| (left) and Mcllwain’s L-shell (right); the black 
cnrves are a £t of points with eqnal percentages of galactic and albedo protons, while the 
red cnrve denotes the Stormer vertical cntoff for the PAMELA epoch. Bottom panels: 
corresponding rigidity prohles, for different valnes of |A| (left) and L (right); valnes at bin 
center are reported in labels. Lines are to gnide the eye. 
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Fig. 7.— Top panel; comparison between fit values obtained with (red) and without 

(blue) the inclusion of external geomagnetic field sources, as a function of L-shell. Bottom 
panel; corresponding relative difference in percentage. 
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Fig. 8.— Mean differential energy spectra (GeV“^ sr“^) in the SAA region {B < 

0.23 G). Resnlts for the different proton populations are shown: quasi-trapped (QT, blue), 
precipitating (UT^, green) and pseudo-trapped (UT^, red). 
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